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ABSTRACT 

The blast-wave model for gamma-ray bursts (GRBs) has been called into question 
by observations of spectra from GRBs that are harder than can be produced through 
optically thin synchrotron emission. If GRBs originate from the collapse of massive 
'-^ | stars, then circumstellar clouds near burst sources will be illuminated by intense 7 

■ radiation, and the electrons in these clouds will be rapidly scattered to energies as 

■ large as several hundred keV. Low-energy photons that subsequently pass through the 
hot plasma will be scattered to higher energies, thus hardening the intrisic spectrum. 
This effect resolves the "line-of-death" objection to the synchrotron shock model. 
Illuminated clouds near GRBs will form relativistic plasmas containing large numbers 

\ of electron-positron pairs that can be detected within ~ 1-2 days of the explosion 

before expanding and dissipating. Localized regions of pair annihilation radiation in 

■ the Galaxy would reveal past GRB explosions. 
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52 ' 1- Introduction 

C3 ■ 

The identification of flaring and fading X-ray, optical and radio counterparts to gamma-ray 
burst (GRB) sources (e.g., Costa et al. 1997; van Paradijs et al. 1997; Djorgovski et al. 1997; Frail 
et al. 1997), and the large energy releases implied by redshift measurements, find a consistent 
explanation in an expanding relativistic blast-wave model (Paczyhski & Rhoads 1993; Meszaros 
& Rees 1997). As a result of Beppo-SAX and optical follow-on observations, the redshifts of 
about one dozen GRBs with durations greater than ~ 1 s have been measured. The distribution 
of redshifts is broad and centered near z ~ 1, corresponding to the cosmological epoch of active 
star formation (Hogg & Fruchter 1999). GRBs are extremely luminous and energetic at hard 
X-ray and 7-ray energies. The degree of GRB collimation is unknown, but peak directional 7-ray 
luminosities and energy releases as large as dL/d£l ~ 3 x 10 51 ergs (s-sr)" 1 and dE/dQ ~ 3 x 10 53 
ergs sr _1 , respectively, have been measured (Kulkarni et al. 1999). Less powerful GRBs and less 
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luminous episodes during the GRB produce smaller 7-ray powers, but the apparent isotropic 
7-ray luminosities from typical GRBs could regularly reach values exceeding 10 50 L5o ergs s" 1 with 
L50 ~ 1, with some GRBs reaching L50 > 10 2 . Because the energy radiated in 7 rays is less than 
the total energy released by a GRB, the apparent isotropic energy release of GRB sources could 
often reach values of 10 54 £ ; 54 ergs, with E54 ~1. 

Cosmological gamma-ray burst and afterglow observations are best explained through the 
fireball/blast- wave model, where the deposition of large quantities of energy into a small region 
yields a fireball that expands until it reaches a relativistic speed determined by the amount of 
baryons mixed into the fireball (see, e.g., Piran 1999 for a review). Nonthermal synchrotron 
radiation from energetic electrons in the relativistic blast wave is thought to account for the 
origin of the prompt 7-ray emission and afterglow radiation. This paradigm has been called into 
question, however, by observations of very hard X-ray emission during the prompt 7-ray luminous 
phase of a significant number of GRBs (Crider et al. 1997; Preece et al. 1998). Photon fluxes 
4>{e) oc e~ ax with ax ~ 0, where e = hu/m e c 2 is the dimensionless photon energy, have been 
observed in 5-10% of GRBs that are bright enough to permit spectral analysis. This strongly 
contradicts the optically-thin synchrotron shock model, which predicts that only radiation spectra 
with ax > 2/3 can emerge from the blast wave. In view of the severity of this challenge to the 
model, these observations have been termed the "line-of-death" to the synchrotron shock model. 
Possible explanations for this phenomeonon involve photoelectric absorption by optically thick cold 
matter (Liang & Kargatis 1994; Brainerd 1994; Bottcher et al. 1999), synchrotron self-absorption 
(Crider & Liang 1999, Granot, Piran, & Sari 2000, Lloyd Sz Petrosian 1999), Compton scattering 
(Liang 1997; Liang et al. 1999), or the existence of a pair-photosphere (Meszaros & Rees 2000) 
within the blast wave. Except for the last model cited, these explanations are inconsistent with 
the standard synchrotron shock model. Here we offer a solution to this problem that is consistent 
with the standard model and recent observations pointing to a massive star origin of GRBs. 

2. Massive Star Origin of GRBs 

Considerable evidence linking the sources of GRBs with star-forming regions has recently 
been obtained (e.g., Lamb 1999). For example, the associated host galaxies have blue colors, 
consistent with galaxy types that are undergoing active star formation. GRB counterparts are 
found within the optical radii and central regions of the host galaxies (e.g. Bloom et al. 1999a), 
rather than far outside the galaxies' disks, as might be expected in a scenario of merging neutron 
stars and black holes (Narayan, Paczyhski, & Piran 1992). Lack of optical counterparts in some 
GRBs could be due to extreme reddening from large quantities of gas and dust in the host galaxy. 
This, together with the appearance of supernova-like emissions in the late time optical decay 
curves of a few GRBs (e.g., Bloom et al. 1999b) and weak X-ray evidence for Fe K Q -line signatures 
(Piro et al. 1999), supports a massive star hypernova/collapsar (Woosley 1993; Paczyhski 1998) 
origin for the long duration gamma-ray bursters. 
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The observations thus favor a model for GRBs involving the collapse of the core of a £ 3OM 
star to a black hole, with the collapse events producing fireballs and relativistic outflows with large 
directed energy releases. Earlier treatments of the blast-wave model considered systems where 
the density of the surrounding medium is either uniform or monotonically decreasing as a result 
of a circumstellar medium formed by a hot stellar wind (Meszaros, Rees, k, Wijers 1998). Until 
recently (Chevalier & Li 1999; Li & Chevalier 1999), less attention has been paid to the actual 
environment found in the vicinity of massive stars. For this we consider rj Carinae (Davidson & 
Humphreys 1997), the best-studied massive star that might correspond to a GRB progenitor. 
It is an evolved star with present-day mass > 90M©, distance of 2300 ± 200 pc, and lifetime of 
about 3 million years. It anisotropically ejects mass at a current rate of ^ 0.003M Q yr _1 to form 
its unusual "homonculus nebula." Several Solar masses of material surround 77 Carinae. In the 
immediate vicinity of the central star, dense clouds of slow-moving gas with radii r ~ 10 15 cm and 
densities between 10 7 and 10 10 cm -3 were discovered with speckle techniques (Hofman & Weigelt 
1988) and confirmed with high-resolution HST observations (Davidson et al. 1995). This material, 
moving with speeds of ~ 50 km s _1 , is apparently ejected nonuniformly from the equatorial zone, 
but may remain trapped by the gravitational field of the star. Inferences (Davidson & Humphreys 
1997) from [Fell] observations suggest that £ O.O2M of gas are contained within ~ 2 x 10 16 cm, 
implying a volume-averaged gas density ^ 7 x 10 5 cm~ 3 . Model results (Bottcher 1999) imply 
that a dense (n ~ 10 12 cm" 3 ) torus of gas at mean distance d ~ 2 x 10 15 cm and with a 10-fold 
enhancement of Fe relative to Solar abundance is required to explain the Fe K a emission weakly 
detected from GRB 970508 with Beppo-SAX (Piro et al. 1999). 

Guided by the optical observations of 77 Carinae, we assume that the volume-averaged density 
of gas at d £ 10 16 cm of a GRB source is (n) = 10 6 ri6 cm~ 3 . Dense clouds of radius 10 15 ri5 cm 
and radial Thomson depths tt = n c arr are assumed to be embedded within this region, so that 
the mean density of particles in a cloud is n c = 1.5 x 10 9 tt/Vi5 cm~ 3 . Thus tt ~ 1 clouds located 
very close to a GRB source are consistent with the observations of dense blobs near r\ Carinae. 
The deceleration length scale of a blast wave with initial Lorentz factor Tq = 100F2 in a uniform 
medium is x d = (3Eq/AttTq (n)mpC 2 ) 1 / 3 = 2.5 x 10 15 (E^/T^uq) 1 ^ cm; hence the blast wave would 
emit a significant fraction of its energy before reaching distances of ~ 10 16 cm. The deceleration 
time, which corresponds to the duration of the prompt 7-ray luminous phase of a GRB in the 
external shock model (Rees & Meszaros 1992), is t d = (1 + z)x d /(cTl) = 8(1 + z)(E 54 /T%n e ) 1/3 
s. These parameters are not unique, and we expect that GRBs display a wide range of energies, 
Lorentz factors and surrounding mean densities that could accommodate the diverse range of 
GRB observations. 

3. Blast- Wave/Cloud Interaction 

A wave of photons impinging on a cloud located 10 16 di6 cm from the explosion center 
will photoionize and Compton-scatter the ambient electrons to energies characteristic of the 
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incident 7 rays (Madau & Thompson 1999). The 7-ray photon front has a width of ~10-100 lt-s, 
corresponding to the duration of the GRB, whereas the plasma cloud has a width of ~ 3 x 10 4 ri5 
lt-s, so that the radiation effects must be treated locally. The radiation force driving the electrons 
outward is balanced by strong electrostatic forces from the more massive protons and ions that 
anchor the system until the net impulse is sufficient to drive the entire plasma cloud outward. The 
Compton back-scattered photons provide targets for successive waves of incident GRB photons 
through 77 pair-production interactions (Thompson & Madau 1999). Higher-energy photons are 
preferentially attenuated, forming an additional injection source of ^ 1 MeV electron-positron 
pairs. The nonthermal electrons and pairs will Compton scatter successive waves of photons, 
thereby modifying the incident spectrum. The pairs, no longer bound by electrostatic attraction 
with the ions, will be driven outward by both radiation forces and restoring electrostatic fields 
to form a mildly relativistic pair wind passing through the more slowly moving normal plasma. 
Shortly after the 7-ray photon front has passed, the decelerating blast wave from the GRB will 
plow into the cloud, shock-heating the relativistic plasma. 

Nonthermal synchrotron photons with energy e impinge on the atoms in the cloud with a flux 
which can be parametrized as 

$(e) = (iird 2 )- 1 £2- [ , - ._.„ 1 - . - ] (1) 

m e c 2 e 2 d L (e/e ) 2 / 3 + (e/e )^ J 

(Dermer, Chiang, & Bottcher 1999), where eo ~ 1 is the photon energy of the peak of the vF v 
spectrum, a 7 ~ 2-3 is the photon spectral index at energies e 3> eo, and £1 = [3/4 + (a 7 — 2) — 1 ]. 
Hydrogen, the most abundant species in the cloud, will be ionized on a time scale of 
5 x 10~ 5 (1 + z)d1 e CieQ^ 3 1 X50 s. Fe features might persist briefly during the early periods of 
very weak GRBs on a time scale of 4 x 1CT 3 (1 + z)d1 6 de^ 3 /L50 s, and would be identified by 
a rapidly evolving Fe absorption feature at 9.1/(1 + z) keV. After the H and Fe are ionized, the 
coupling between the GRB photons and gas is dominated by Compton scattering interactions. 
Pair production through photon-particle processes are negligible by comparison with Compton 
interactions except for photons with e it 200. A lower limit to the time scale for an electron to be 
scattered by a photon is tr(s) ~ 15(1 + z)(if 6 Cieo/C2-^50) where C2 = [3 + (a 7 — l) -1 ], assuming 
that all Compton scattering events occur in the Thomson limit. The Klein-Nishina decline in the 
Compton cross section will increase this estimate by a factor of ~ 1-3, depending on the incident 
spectrum. Most of the electrons in the cloud will therefore be scattered to high energies during a 
very luminous (L50 S> 1) GRB, or when the cloud is located at d\Q <C 1. 

The average energy transferred to an electron at rest when Compton-scattered by a photon 
with energy e is Ae = e 2 /(l + 1.5e) (this expression is accurate to better than 18% for e < 10 3 ). 
Defining rj = 7 — 1 as the dimensionless electron kinetic energy, we can easily estimate the 
production rate f{rj) of electrons scattered to energy rj in the nonrelativistic (rj <C 1) and 
extreme relativistic (rj 3> 1) limits, noting that f(rj)dr) oc <E>(e)<7c(e)cfe and letting rj = Ae. In 
the former limit, the Compton cross section crc(e) — ► ax and 3>(e) oc e~ 2 / 3 , so that f(rj) oc ry -5 / 6 
when rj <C min(l, eg). In the high energy limit, C7c(e) oc ln(3.3e)/e and $>(e) oc e"" 7 , so that 
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f(rj) oc ln(2.2r/) when r\ ^> max(l,eg). Thus electrons are Compton-scattered on 

the time scale derived above to form a hard spectrum that turns over at kinetic energies of 
^ 500xmin(l, eg) keV. For a GRB with eo ~ 1, most of the kinetic energy is therefore carried by 
nonthermal electrons with energies of ~ 500 keV. Successive waves of photons that pass through 
this plasma will continue to Compton-scatter the nonthermal electrons. Only the lowest energy 
photons, however, will be strongly affected by the radiative transfer because both the Compton 
scattering cross section and energy change per scattering is largest for the lowest energy photons. 

Following the initial wave of photons, successive photon fronts also encounter the back- 
scattered radiation (Madau &; Thompson 1999; Thompson & Madau 1999). The kinematics of the 
Compton process dictate that the energy e s of a photon back-scattered through 180° by an electron 
at rest is e s = e/(l + 2e); thus e s cannot exceed 1/2 the electron rest-mass energy. Head-on collisions 
of the back-scattered photons by primary GRB photons with e x > 2/e s = 2 + 2^3 can thus 
produce nonthermal e + -e~ pairs. The cross section for 77 pair production peaks near threshold 
with a value of ~ 0"t/3. The 77 pair-production optical depth T 77 (ei) of a photon that trails the 
onset of the GRB by At seconds can be estimated by noting that the photon traverses a distance 
~ r through a backscattered radiation field of spectral density n s {e s ) » n e OT ■ Ai • <5[e s /(l — 2e s )] 
- a more accurate calculation would replace the term At • $[e s /(l — 2e s )] by an integral over 
the time-varying flux. Approximating r 77 (ei) r(or/3)Ae s n s (2/ei), where Ae s ~ 2/ei is the 
bandwidth that is effective for producing pairs, we obtain 

_ fc r rAt[s] L 50 fc(ei) 1 

TTr(ei) ~ °-° 2 d^Ci [ VJ^TVJ^ ] ' (2) 
where e' = 2/(ei — 4). The coefficient results from a more detailed derivation, and the term 
k(e\) = 1 — 4e 1 ' 1 + ei/(ei — 4) is a Klein-Nishina correction. Eq. @) shows that photons with 
energies above several MeV will be severely attenuated in Thomson thick clouds if L50 3> 1 or 

— 1/3 

d\Q <C 1. Photons with MeV energies are most severely attenuated, and r 77 (ei) oc e 1 at energies 
ei S> max(l,eo). The 77 pair injection process provides another source of nonthermal leptons 
with 77 ~ 1. The pairs will not, however, be electrostatically bound but will be accelerated by the 
photon pressure and electrostatic field. 

Fig. 1 shows Monte Carlo simulations of radiation spectra described by eq. (|l]) that pass 
through a hot electron scattering medium. For simplicity, we approximate the hard nonthermal 
electron spectrum by a thermal distribution with temperatures of 100 and 300 keV. These 
calculations show that the lowest energy photons of the primary synchrotron spectrum are most 
strongly scattered, and that the "line-of-death" problem of the synchrotron shock model of GRBs 
(Preece et al. 1998) can be solved by radiation transfer effects through a hot scattering cloud with 
t t Z 1-2. 

According to this interpretation, GRBs displaying very hard spectra could display one break 
from the intrinsic synchrotron shock spectrum and a second break from the scattering process. In 
the examples shown in Fig. 1, these two breaks are so close to each other that they appear as one 
smooth turnover. Two breaks are observed from GRB 970111 (Crider & Liang 1999), a GRB that 
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strongly violates the "line of death." A prediction of this model is that GRBs showing such flat 
X-ray spectra should also display softer MeV spectra than typical GRBs due to 77 attenuation 
processes in the hot scattering cloud. 

4. Observational Signatures of the Flash-Heated Cloud 

The Compton-scattered electrons transfer momentum to the N p = 4 x 10 54 r^ 5 ng protons 
of the cloud through their electrostatic coupling. If the radiation efficiency is £ r , then the 
Compton impulse gives each proton in the cloud ps £ r [l — exp(— TT^Eitr 2 / (And 2 N p ) « 
40 [1 — exp(— T'r)](^r/0.1)-E ; 54/(df 6 ri5n9) MeV of directed energy. Pairs, by contrast, will be 
accelerated to mildly relativistic speeds until Compton drag or streaming instabilities limit further 
acceleration. In the simplification that the medium interior to the cloud is uniform, and neglecting 
pair- loading of the swept-up material (Thompson & Madau 1999), the decelerating blast wave 
follows the dynamical equation T(x) = T${x/ 'xd)~ 9 ', where g = 3/2 and 3 for adiabatic and 
radiative blast waves, respectively. Using the standard parameters adopted here, the blast wave 
slows to between £ = 0.01 — 0.1 of its initial Lorentz factor before reaching a cloud at d = 10 16 
cm. Even considering the radiative acceleration of the cloud, the blast wave reaches the cloud 
at time tb w = td{d/xd)^ 29+1 ^ / (2g + 1) rS idyn, where the dynamical time scale of the cloud is 
^dyn = t/c = 3 x 10 4 ri5 s. Because the cloud is so dense, a large fraction of the residual energy 
of the blast wave is deposited into the N p particles of the cloud. Thus each proton in the cloud 
receives an additional m p (3 2 c 2 « £Eirr 2 / (And 2 N p ) « 40(f/0. 1)E 54 /(d 2 l6 r 15 n 9 ) MeV of kinetic 
energy, divided roughly equally into directed outflow and random thermal energy. 

If the circumstellar medium at d\Q » 1 is much more dilute than the interior region, 
as suggested by observations of r] Carinae, then we can neglect further interactions of the 
cloud/blast wave system with their surroundings. The observational signatures and fate of the 
cloud at late times can be outlined by comparing time scales. The cloud expands on a time scale 
^cx = tdyn/ = 1-5 x 10 5 ri5/(/3/0.2) s. The basic time scale governing radiative processes in the 
cloud is the Thomson time <t = (n e <7Tc) _1 = 5 x 10 4 /ng s = idyn/ r T- The electrons thermalize on 
a time scale tx/lnA <C tdyn, where the Coulomb logarithm In A « 20. The protons transfer their 
energy to the electrons on the time scale t ep = 3 l 2 (m p /m e )t^ / hi A, where = kT/m e c 2 is an 
effective dimensionless electron temperature, and we assume collective plasma processes for energy 
exchange to be negligible. 

The flash-heated cloud can evolve in two limiting regimes. When the external soft-photon 
energy density is small, the system emits by a hard bremsstrahlung spectrum with ~ 1 
and luminosity Lff = N e afm e c 2 l l 2 /ti ~ 5 x 10 41 rf 5 ng# 1//2 ergs s _1 . In the more likely 
case when abundant soft photons are present, for example, from the Compton echo (Madau, 
Blandford, & Rees 1999), then Compton cooling will balance ion heating to produce a 
luminous Comptonized soft-photon spectrum with effective temperature ~ 0.1 and luminosity 
L c = iEirr 2 /{And 2 t ep ) w 2 x 10 45 (£/0.1)S 5 4rf 5 n 9 d 1 " 6 2 (6'/0.1)- 3 / 2 ergs s -1 . In either case, the 
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spectra persist until the plasma expands and adiabatically cools, that is, for a period ~ t ex ~ day. 
The hot bremsstrahlung plasma will be too dim to be detectable with current instrumentation, but 
a 50-100 keV Comptonized plasma at redshift z ~ 1 and luminosity distance of 10 28 I?28 cm would 
have a flux of ~ 1.6 x 10" 12 (C/0.1)S 5 4rf 5 n 9 (i^ 6 2 ((9/0.1)- 3 / 2 L' 28 2 ergs cm~ 2 s -1 . The hot plasma 
formed by a nearby GRB at z ~ 0.1 would be easily detectable with the INTEGRAL and Swift 
missions at hard X-ray and soft 7-ray energies. In either case, e + -e~ pairs would be formed with 
moderate efficiency, and the cooling, expanding plasma would produce a broad pair annihilation 
feature (Guilbert & Stepney 1985). The residual pairs formed in the relativistic plasma and the 
pair wind would diffuse into the dilute interstellar medium with density nisM to annihilate on a 
time scale (^isMfTc)" 1 ~ 2 x 10 6 /nisM yr. If the energy intercepted by a single cloud is converted 
to pairs with a conservative 1% pair yield, past GRBs in the Milky Way would be revealed by 
localized regions of annihilation radiation with flux ~2x 10- 5 £ 54 nisM (d/lOkpc)" 2 0.511 MeV ph 
cm -2 s _1 . The high-latitude annihilation feature discovered with OSSE on the Compton Gamma 
Ray Observatory (Purcell et al. 1997), or other localized hot spots of annihilation radiation that 
will be mapped in detail with INTEGRAL, could reveal sites of past GRB explosions. 

CD thanks B. Paczyhski for stressing the importance of massive star observations in 
developing blast-wave models of GRBs. The work of CD is supported by the Office of Naval 
Research and NASA Astrophysical Theory Program (DPR S-13756G). The work of MB is 
supported by NASA through Chandra Postdoctoral Fellowship grant PF 9-10007, awarded by 
the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for 
NASA under contract NAS 8-39073. 
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Fig. 1. — Radiation transfer effects on GRB emission that passes through electrons energized by 
an earlier portion of the photon front. The intrinsic spectrum eq.(l), with ax = 2/3, a 7 = 2.5 and 
eo = 0.5, is shown by the thick dashed curve. The nonthermal electrons are approximated by a 
thermal distribution with temperatures of 100 keV (thin curves) and 300 keV (thick curves), and 
Thomson depths tt = 1 (solid curves) and tt = 2 (dotted curves). Spectral indices ax calculated 
at 50 keV are 0.5 (T = 100 keV, r T = 1), 0.44 (T = 300 keV, r T = 1), 0.14 (T = 100 keV, r T = 2), 
and 0.05 (T = 300 keV, t t = 2). 



